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We consider general relativity with a cosmological constant as a perturbative expansion around a 
completely solvable diffeomorphism invariant field theory. This theory is the A oo limit of general 
relativity. This allows an explicit perturbative computational setup in which the quantum states of 
the theory and the classical observables can be explicitly computed. The zeroth order corresponds 
to highly degenerate space-times with vanishing volume. Perturbations give rise to space-times with 
non-vanishing volumes in a natural way. The spectrum of area- and volume-related observables 
constructed by coupling the theory to matter can be directly assessed. An unexpected relationship 
arises at a quantum level between the discrete spectrum of the volume operator and the allowed 
values of the cosmological constant. 



The recent years have been very fertile for the de- 
velopment of canonical quantizations of general relativ- 
ity. There are now available more than one highly 
non-trivial mathematically consistent, finite, anomaly- 
free theories of canonical quantum gravity. Instrumental 
in these developments have been the underlying advances 
in mathematical techniques for dealing with infinite di- 
mensional nonlinear spaces, like the theory of cylindrical 
functions and associated measures and the introduc- 
tion of spin networks to eliminate the over-completeness 
of the Wilson loop basis Q (for a recent summary, see 
the review article by RovcUi [||). In spite of the ad- 
vances, most of the results obtained from these theories 
up to now concern statements made at a kinematical level 
(without imposing the Hamiltonian constraint). One of 
the remaining main challenges is to obtain results that 
hold at a dynamical level, imposing the full set of con- 
straints of the theory. This requires making the theories 
computable, in the sense of obtaining classical observ- 
ables to promote to quantum observables and solving for 
the physical set of states that are annihilated by all the 
constraints. Only in this setting will one be able to make 
contact with trustworthy physical predictions. 

The fact that implementing the full dynamics of canon- 
ical quantum gravity is hard to do has plagued the sub- 
ject since its inception in the 60's by DeWitt For in- 
stance, one does not know even at a classical level how to 
compute observable quantities for canonical general rel- 
ativity (quantities that are invariants under the symme- 
tries of the theory or — in the canonical language — that 
commute with the constraints), and it is known that such 
quantities are potentially quite involved Q . As a conse- 
quence, it is very hard to make contact with the semi- 
classical limit, since one does not have at hand observable 
quantities to compute. Moreover, the space of quantum 
states that are annihilated by the quantum constraints 
cannot be straightforwardly computed . None of these 
issues was directly tackled by the significant amount of 



progress that was made in terms of the introduction of 
connection type variables and loop representations for 
quantum gravity that we mentioned at the beginning. 
There is general consensus that an approximation scheme 
is desirable in such a way as to make the above problems 
solvable. The challenge consists in finding an approxima- 
tion scheme that will not conflict with non-perturbative 
nature of the canonical quantization, in particular the 
diffeomorphism invariance of the theory. 

In this paper we would like to propose a computational 
scheme that allows to explicitly address the issue of find- 
ing observables for the theory (at a classical level) and 
to find quantum states annihilated by all the constraints. 
The central idea will be to consider general relativity as 
a power series expansion in terms of the inverse cosmo- 
logical constant, and to view the full theory as a "per- 
turbation" of a zeroth order solvable theory. The zeroth 
order theory corresponds to general relativity in the limit 
in which the cosmological constant A —^ oo. The use of 
an approximation based on a large cosmological constant 
has been pioneered in a different context (that of study- 
ing the spectrum of the physical Hamiltonian obtained by 
"de-parameterizing" the theory) by Smolin and also 
by Rovelli [|lO|, with some level of success. We will how- 
ever be able to make further progress, in significant part 
due to the new maturity of the field we mentioned above 
but also because we push the perturbative approach at all 
levels: constraints, observables and equations of motion, 
both classically and quantum mechanically. It might be 
argued that a scheme where the cosmological constant 
is large (in Planck units) lacks physical relevance, given 
is small present value. However, having a regime of the 
theory in which it is completely solvable can be a valu- 
able tool, even if such regime is unphysical. Moreover, 
as we will discuss below, there is some evidence that our 
approach might be useful beyond its apparent domain of 
validity. It might be the case that the situation is analo- 
gous to the large N approximation in QCD or the strong 
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coupling expansions of lattice field theory, that a priori 
look less promising that what they end up being. 

The starting point is to consider canonical quantum 
gravity with a cosmological constant. The Hamiltonian 
constraint of the theory is. 



H{N) = Ho{N) + A / d^xN{x)^/detq{x) 



(1) 



where N{x) is a smearing function, Ho{N) is the Hamil- 
tonian constraint without a cosmological constant and 
detg is the determinant of the spatial metric. If one now 
considers the limit A — ^ oo, and re-scales the constraint 
by 1/A [0 one is left with a theory, which we will call 
"zeroth order" theory for which the Hamiltonian con- 
straint is, ij(o)(iV) = / d^xN{x)y^detqix), that is, the 
Hamiltonian constraint is just the square root of the de- 
terminant of the spatial metric. In addition to this the 
theory has the ordinary diffeomorphism constraint and 
if one uses Ashtekar variables there will also be a Gauss 
law constraint (both are independent of A). Imposing 
classically the Hamiltonian constraint of the zeroth or- 
der theory, one immediately sees that it corresponds to 
metrics of identically vanishing determinant, that is, de- 
generate metrics. 

The dynamics of the zeroth order theory is easily 
solved. One can readily construct quantities that have 
vanishing Poisson bracket with the Hamiltonian con- 
straint of the zeroth order theory. For instance, one can 
consider any function depending only on the three metric 
(or if one is using Ashtekar variables (i?f , A^), the densi- 
tized triads Ef) and not on its canonically conjugate mo- 
menta. To be a genuine observable the quantity should 
also have vanishing Poisson bracket with the diffeomor- 
phism constraint. A way of constructing quantities with 
this property is to couple the theory with matter fields 
and physically define geometric objects like surfaces or 
volumetric regions through properties of the matter fields 
and then compute the areas or volumes of these 
regions. The usual problem with this construction is that 
the addition of matter adds terms to the Hamiltonian 
constraint as well, and the areas and volumes computed 
as we just described do not commute with the Hamilto- 
nian constraint. However, in the A — > co limit the matter 
terms drop out of the Hamiltonian and one ends up (in 
the zeroth order theory) with genuine observables (this 
is similar to what happens in models without a Hamilto- 
nian constraint as was discussed by Rovelli in the Husain- 
Kuchai^ model ||l^). A concrete example of an ob- 
servable for the zeroth order theory is obtained by cou- 
pling the theory to an antisymmetric tensor field with the 
standard Kalb-Ramond free Hamiltonian as proposed by 
Smolin jl^ . The antisymmetric tensor field is naturally 
associated with surfaces of which one can compute the 
area. The resulting quantity, which essentially is a con- 
traction of the triads with the canonically conjugate mo- 
menta of of the Kalb-Ramond field, is naturally diffeo- 



morphism invariant. One can build a "loop-surface" rep- 
resentation (see 14 for details) for the coupled theory 
by expanding the gauge invariant wavefunctions in a ba- 
sis of states |E, s > labeled by equivalence classes under 
diffeomorphisms of spin networks s and three-surfaces E. 
A natural generalization of the Ashtekar-Lewandowski 
Q measure allows to introduce an inner product. On 
this basis the area defined by the fields has a discrete 
spectrum. 

The Hamiltonian constraint of the zeroth order theory 
can be readily realized in the kinematical space |^,|l^ of 
non-diffeomorphism invariant wavefunctions as, 

H^"\N)\T,,s>^ liiJiy^ N{v)V{v,e)\J:,s > (2) 

where we have represented, following Thiemann ||l[, the 
determinant of the metric as the volume operator of an 
infinitesimal region of size e centered on the vertices v 
of the spin network. The operator is independent on 
e on this space of states and therefore the above limit 
is well defined. To consider the action on diffeomor- 
phism invariant states one evaluates the inner product 
< {S'},{s}|i?(°)(A^)|S,s > which is well defined. We 
can readily find the space of states that are annihilated 
by the Hamiltonian constraint of the zeroth order the- 
ory: they correspond to states based on spin networks 
such that the localized volume operator vanishes on each 
of their vertices. For trivalent vertices this is the case 
for all states. For four and higher valent vertices the op- 
erator vanishes automatically if the vertices are planar 
(this is not true for the volume operator of [Q, but it 
holds for the definition of ||l^). If the vertices are not 
planar it is also possible to find values of the valences of 
the incoming edges such that the operator vanishes. We 
collectively denote all these states as |S, s° >. Since the 
space of solutions of the Hamiltonian constraint is a sub- 
space of the states we were considering, it can be readily 
endowed with an inner product which is just the restric- 
tion to this space of the inner product derived from the 
Ashtekar-Lewandowski measure. Notice that the area 
observables we introduced are well defined on this space. 

We have therefore seen that the zeroth order theory is 
well defined. It is non-perturbatively quantized in a basis 
of spin networks and surfaces, one can solve for the phys- 
ical states, construct observables, and introduce a physi- 
cal inner product. This is therefore a sound starting point 
for a perturbative expansion. The resulting perturbative 
theory is remarkably similar to bound state perturbation 
theory in quantum mechanics in the degenerate case. The 
operator 7?^°^ (which can be readily diagonalized in the 
basis of states considered) has a degenerate spectrum. 

To consider the higher order corrections to the zeroth 
order theory, one is interested in solving the eigenvalue 
problem 



< 0(A-i)|ij("' + A-^H^^'I = e(A) < 0(A- 



(3) 



2 



with eigenvalue e(A) equal to zero up to the desired order 
in perturbation theory. We will assume that the states 
and the eigenvalue can be expanded in power series in 
the inverse cosmological constant, and address the re- 
sulting equations order by order. Since the spectrum of 
the zeroth order Hamiltonian is degenerate, we apply de- 
generate perturbation theory to solve this problem. One 
starts by considering a state with eigenvalue e^^"^ for the 
zeroth order Hamiltonian constraint, that is, a spin net- 
work state with a definite volume. Notice that this is not 
a state of the zeroth order theory. Considering first order 
corrections and solving the equations resulting from (|^) 
one determines the first order corrections and a polyno- 
mial e(A~^) which, set to zero, determines the value of 
A for which the perturbative treatment is valid. That is, 
given an eigenvalue of the zeroth order Hamiltonian, one 
constructs a solution of the first order theory that has 
vanishing Hamiltonian for a particular value of lambda. 
That is, as a consequence of the discrete spectrum of the 
zeroth order Hamiltonian inherited from the discreteness 
of the volume operator, the cosmological constant is, in 
this approach, quantized. 

This might sound unusual, but imposing the constraint 
perturbatively is analogous to considering an atom in a 
magnetic field and trying to find the states with a given 
energy level. Such energy level can be reached by choos- 
ing a nearby energy state of an unperturbed atom and 
then tuning the magnetic field to a given value. Various 
states could correspond to the same energy, depending on 
the value of the magnetic field. The approximate value of 
the magnetic field required will also depend on how high 
one goes in perturbation theory. A detailed discussion of 
how this is implemented in several systems can be seen 
in 0. 

A remarkable aspect is that all the calculations in- 
volved in the perturbative approach are completely ex- 
plicit when one works in the spin network basis. One 
starts by considering eigenvalues and eigenvectors of the 
zeroth order Hamiltonian. Since this Hamiltonian is the 
volume operator, which is well understood, the space of 
states is under control (even when it is non- vanishing) . 
The equation that determines the first order correction 
involves evaluating the Hamiltonian constraint on the 
states with a given volume. Such calculations are com- 
pletely explicit given, for instance, the Hamiltonian in- 
troduced by Thiemann or the Hamiltonian in the 
Vassiliev arena . The only subtle element in the calcu- 
lation is that given the infinite degeneracy of the volume 
operator, the first order correction might involve infinite 
superpositions and therefore may lead us to corrections to 
the kinematical inner product. It is worthwhile noticing 
that in general the first order corrections have non-zero 
volume, even if one started from the zeroth order state 
with vanishing volume, as long as one did not restrict to 
either planar or trivalent spin networks. So the resulting 
theory is a theory of non- vanishing volume space-times. 



The above construction has a classical counterpart 
when one tries to find observables for the theory. An 
observable for the first order theory can be written as a 
perturbation of the observables of the zeroth order level 
theory we discussed above, 

Oa(^,^) -0^°Hvr,^)+A-iO(i)(^,^). (4) 

One would like to request that these observables have 
vanishing Poisson brackets with the Hamiltonian of the 
theory (more precisely, weakly vanishing, but the exten- 
sion to such case is immediate) . Expanding such require- 
ment in powers of A~^ one gets to zeroth and first order 
in A-\ 

|o("),i/(0)} ^0 (5) 
{oW,i/("'} + {o("),i/(i)} =0 (6) 

The first equation determines 0*^°-' and the second one 
leads to a (functional) linear partial differential equation 
for O^^). The construction can be readily continued to 
higher orders. In all cases one obtains a linear partial dif- 
ferential equation, albeit with a more and more complex 
inhomogeneous term. It should be noted that one can 
obtain many observables starting with different solutions 
to the first equation. The linear partial differential equa- 
tions are not hard to solve, given that the coefficients of 
the derivatives are functions of the triads, whereas the 
derivatives are with respect to the connections. Given 
its simplicity, the system is always integrable and there- 
fore yields all the observables of the theory. The solu- 
tions to these equations will generically be ill defined or 
non-convergent in certain regions of phase space. This 
is expected in theories with complicated dynamics and is 
not a pathology of the method per se but rather a feature 
of the observable quantities. In particular, one might be 
concerned that the current approach produces too many 
observables. For instance, there exist pathological sys- 
tems that have less observables than degrees of freedom 
(e.g. the heavy asymmetric top). A priori our method 
seems to yield in these cases extra observables, but on 
closer examination the quantities constructed are not ac- 
ceptable (in the case of the top they take values on a 
restricted portion of phase space not preserved by evolu- 
tion). We discuss several pathological examples and how 
the method correctly handles them in . 

Let us consider, as an simple example of interest in 
gravity, the computation of an observable for Bianchi 
cosmologies. Observables for these models (with cosmo- 
logical constant) had never been studied before. We will 
write the Bianchi models in terms of Ashtekar's variables 
following the notation of , in which one codes all the 
information in three variables associated with the time 
dependent portion of the triads E^,E'^, and their con- 
jugate momenta Ai^ A2, A3. The only constraint of the 
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theory given the homogeneity and diagonahty of the met- 
ric is the Hamiltonian, which rescaled by the cosmological 
constant reads, 

H = E^E^E^ + A-^ {E^E^{AiA2 - eA^) + cycHc) , (7) 

where e = 1 corresponds to the Bianchi IX model and 
e = to the Bianchi I model and the extra terms are 
obtained by cyclically permuting the indices 1,2,3. It is 
straightforward to repeat the construction for the other 
class-A Bianchi models. 

The observables of the zeroth order theory are given 
by functions 0° = F{E^,E^,E^,E'^Ai - E^A2,E^Ai - 
E^As) (if f one considers FiE'^Ai - E^A2,E^Ai - E^A^) 
these objects are exact observables to all orders for the 
Bianchi I case). 

As an example, we can start with an observable of the 
zeroth order theory O*^"^ = E^. The condition (||) trans- 
lates into the following partial differential equation 



dAi 



■E^E^ - E^E^Ao 



E^E^A:, 



eE^E- 
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cyclic = 0. 
(8) 



This yields a first order correction, which modulo the 
Hamiltonian constraint can be written as, 



(1) 



Ai 
' E^E^ 



[-E^E^A2 + {E^)^Ai-E^E^A3), (9) 



up to an arbitrary additive function F{E^ ,E^,E^, E^Ai — 
E^A2,E^Ai - E^A3) that solves the homogeneous part 
of the equation. Expressions for higher orders can eas- 
ily be generated using computer algebra, we do not list 
them here for reasons of space. The behavior of the 
approximate observables can be studied explicitly us- 
ing the exact solution [|8| of the equations of motion 
(in the Bianchi I case). One observes that, evaluated 
on the trajectories of the solution, the coefficients in 
the power series expansion lose their dependence on the 
cosmological constant as an expansion parameter. Care- 
fully choosing the arbitrary functions F one can define 
observables that are finite at the Big Bang and remain 
approximately constant for a significant portion of the 
lifetime of the (recoUapsing) universe. The rate of change 
decreases if one adds further terms in the expansion, ex- 
hibiting convergence, although all information about 
A has disappeared. Therefore quite unexpectedly, the 
method produces solutions that are accurate irrespective 
of the value of the cosmological constant. 

Summarizing, casting ordinary general relativity as a 
perturbative theory in A~^ and starting with background 
theory in which A — > oo allows to explicitly compute 
non-degenerate physical quantum states and observables. 
The technique is applicable in the theory in general and 



also in minisuperspace models. In the general theory, the 
calculations are well defined and tractable, but quite in- 
volved. This opens up a program for the quantization of 
the gravitational field that is well defined computation- 
ally and allows contact with the physical observables of 
the theory, a key ingredient in the road to finding a the- 
ory with a semiclassical regime that makes contact with 
ordinary low energy physics. 
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FIG. 1. The perturbative observable evaluated for a Bianchi I cosmology with negative cosmological constant. The time 
runs from zero, the Big Bang, to tt, the Big Crunch. The four curves correspond to better and better approximations. 
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